Coulomb Drag in Systems with Tunneling Bridges 



CD 

m 



> 

m 
o 

O 



X3 
O 

o 



X 



Yuval Oreg°, and Alex Kamenev'' 
''Department of Condensed Matter Physics The Weizmann Institute of Science 76100 Rehovot, Israel 
''Institute for Theoretical Physics University of California Santa Barbara, CA 93106-4030, U.S.A. 

(February 7, 2008) 

We study the Coulomb drag effect in double layer electronic systems with local tunneling links. 
The possibility of tunneling between the layers leads to a pronounced exchange contribution to 
the transconductance, which is negative and non-vanishing at zero temperature. The diffusive 
motion of the electrons inside each layer in interplay with the electron-electron interaction cause a 
distinguishable singular temperature dependence of the transconductance at low temperatures. 

PADS numbers: 73.20.Dx, 73.40. Gk, 73.50.Yg 



Coulomb drag in spatially separated double-layer elec- 
tronic systems has received a great deal of attention re- 
cently, both experimentally and theoretically [|6|-p^ . 
A setup for studying the drag effect consists of applying 
current to one of the two layers, and measuring the in- 
duced voltage in the second layer. The ratio between the 
two, known also as the transresistance, po, is the main 
characteristic of the drag effect. Theoretically, it is more 
convenient to study the transconductance, ao, which is 
related to transresistance as (for ao ^ cr) 
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where <j = e^Dv is the conductance of a single isolated 
layer (hereafter we discuss the case of identical layers); 
D is the diffusion constant, and v is the density of states 
per unit area. Most of the previous studies were con- 
centrated on systems, in which the coupling between the 
layers is only due to the interlayer Coulomb interaction. 
In this case density fluctuations of the electron liquid 
in the first layer induce fluctuations in the second one, 
which in turn lead to an induced current |]^-[To[|. Since 
it involves classical thermal fluctuations, interacting in 
a Hartree way, the drag effect vanishes at zero temper- 
ature. Moreover, the sign of the effect does not depend 
on the nature of the interaction (repulsion or attraction) 
and has to do only with the sign of the charge carriers 
in the two layers. For example, for carriers of the same 
type, the drag transconductance is strictly positive, i.e., 
the induced current flows in the direction of the applied 
one. 

In this paper we address a qualitatively different mech- 
anism which leads to a transconductance. It takes place 
in the presence of pointlike shortages (or bridges) — 
points where electrons may tunnel between the two lay- 
ers. Such bridges are often present in metallic double- 
layer systems. As we show below the interplay be- 
tween tunneling and the Coulomb interactions lead to 
a transconductance which is non- vanishing at zero tem- 
perature and negative (for carriers of the same charge). 
That is, the induced current flows in the direction op- 
posite to the driving one. The origin of the effect is in 



the exchange interaction between the electrons of differ- 
ent layers, which is possible due to the wavefunctions' 
overlap at the bridges. To illustrate the physics of the ef- 
fect consider a wavepacket propagating in the first layer 
in the direction of an applied current. Once it reaches a 
bridge, a part of the wavepacket tunnels into the second 
layer, while the remaining part continues to move in the 
initial direction. Without interlayer electron-electron (e- 
e) interaction, the wavefront propagating in the second 
layer is spherically symmetric, thus no net current is in- 
duced. (Indeed, after tunneling through a point contact 
the electron completely "forgets" the initial direction of 
its momentum.) In the presence of interaction the wave- 
front moves in a direction that minimized its interaction 
with the initial wavepacket. Thus, for a repulsive inter- 
action the electrons that tunneled move in a direction op- 
posite to that before tunneling. This leads to a negative 
transconductance. As we shall see below such a lowest 
order (in tunneling) mechanism is dominant only at not 
too small temperatures. At lower temperatures the lead- 
ing mechanism involves coherent tunneling to the second 
layer and back to the first one accompanied by intralayer 
Coulomb interactions. The exchange contribution, unlike 
the Hartree term, is not proportional to a small electron- 
hole asymmetry factor ||6|-[l0t|, therefore, the former's ab- 
solute value may well overcome the latter's even for a 
small tunneling rate. In the case of diffusively propa- 
gating wavepackets the exchange transconductance has 
a peculiar temperature dependence, which is described 
below. 

The tunneling rate between the two layers may be de- 
scribed by a mean intralayer lifetime, ri2, which is re- 
lated to the interlayer tunneling conductance per unit 
area, a±, as 



(2) 



We shall assume that the tunneling is weak, i.e., T12 ^ r, 
where t is the mean elastic scattering time within each 
layer. Three energy scales determine the temperature de- 
pendence of the transconductance: h/Ti2 < {Kd)h/Ti2 < 
Ti/t |13[ , where k = 2'Ke^i' is the inverse Thomas-Fermi 
screening radius, and d is the distance between the layers. 
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At high temperatures, Ti/t < T, the motion of the 
wavepackets is balhstic and the exchange contribution to 
the transconductance is negative and temperature inde- 
pendent [its value is discussed below, see Eq. (|l|)]. At 
T < h/r the diffusive character of the electron motion 
should be taken into account. Up to the lowest order in 
the tunneling amplitude we find the following tempera- 
ture dependent contribution to the transconductance 
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The appearance of a contribution which has a singular 
temperature dependence is not entirely unexpected for a 
diffusive system [Q. What is less common is the fact 
that the divergence is so pronounced, ~ , instead of 
being logarithmic, as in the case of Altshuler-Aronov in- 
teraction corrections to the conductivity of 2D systems 
[0. The temperature dependence becomes even more 
singular at smaller temperatures, T < {Kd)h/Ti2 
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The temperature divergences in Eqs. ( |3a| ) and ( |3b| ) are 
cut off by the system size, L, at the Thouless tempera- 
ture, hD/L^. For systems larger than the mean propaga- 
tion distance within one layer, V Dti2 , there is an addi- 
tional temperature range hD/L^ < T < Ti/ti2- Here the 
exchange contribution to the transconductance is dom- 
inated by multiple tunneling processes, leading to the 
form 
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This expression differs from the Altshuler-Aronov inter- 
action correction to the 2D conductivity only by a 
factor of 1/4, which reflects the peculiarity of the drag 
measurement setup. Where the current flows only in one 
half of the system and the potential is measured in the 
other half. 

We now turn to the derivation of Eqs. (|3a|)-(|3c|). The 
Hamiltonian of the system is given by 



H = Hi+H2+H; 
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where ^^1(2) is the Hamiltonian of the first (second) iso- 
lated layer, including elastic disorder, and iJjnt includes 
intcrlayer as well as intralayer Coulomb interactions. The 
first three terms on the r.h.s. of Eq. (^ are traditionally 
involved in the description of the drag effect |p|-p^. We 
add the term describing pointlike tunneling processes 
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where a(a^) and b{b'') are the annihilation (creation) op- 
erators of electrons in the first and the second layers re- 
spectively, and Ti, i = 1 . . . are the random positions 



of N bridges. The tunneling amplitude V is related to 
the lifetime T12 as 2TrvL^N\V\^ = h/ru (cf. Eq. (|)). 




FIG. |l] Two diagrams contributing to the transconductance 
that are second order in tunneling (denoted by x). Full lines 
with arrows are electron Green functions, dashed lines rep- 
resent difFusons, and wavy lines screened interactions. Two 
additional diagrams with arrows in the opposite direction 
should be included. The numbers indicate the layer index. 

At various temperatures the exchange contribution to 
the transconductance originates from the different pro- 
cesses. For example, for {Kd)h/Ti2 < T < Ti/t the dom- 
inant contribution is second order in tunneling and is 
given by the four diffuson diagrams depicted in Fig. ^. 
Performing an analytical continuation and taking the dc 
limit in a standard manner fill] , one finds for the ex- 
change contribution to the transconductance 
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A straightforward summation over the fast electronic mo- 
menta leads to Fxx{uj) — Fxx y^) with 

A- ^ \^ nn^ ^i2(Q,^) 

where Ui2{Q,uj) is the dynamically screened interlayer 
Coulomb interaction. In the diffusive limit {lo < Ti/t) it 
was found in Refs. ||,^ (see also Eqs. (|l^) and (|l^) be- 
low) . Substituting Eq. (^ into Eq. (^ and carefully per- 
forming frequency and slow momenta integrations, one 
recovers Eq. (|3a|). 

At somewhat smaller temperatures the fourth order (in 
the tunneling amplitude) processes may become domi- 
nant. Although being smaller by the additional factor 
of |yp cx 1/ti2, it contains intra- (instead of inter-) 
layer interaction, which is stronger by a Kd > 1 fac- 
tor. Moreover, the fourth order term has a more sin- 
gular temperature dependence (because additional diffu- 
sons poles are involved) which makes it eventually dom- 
inant at T < {Kd)h/Ti2. Examples of the diagrams, con- 
tributing to the fourth order are depicted in Fig. ^. The 
calculation of these diagrams follows the steps outlined 
above and leads to the result given in Eq. (3b). Instead of 
reproducing these calculations we develop a general for- 
malism suitable for multi-tunneling processes, and obtain 
the fourth order contribution as a particular case. 
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FIG. ^ Examples of diagrams contributing to the transcon- 
ductance that are fourth order in tunneling. Diagrams with 
interaction lines connecting "upper" and "lower" Green func- 
tions, as well as diagrams with an opposite direction of elec- 
tron lines should also be included. The numbers indicate the 
layer index. 

The treatment of multiple tunneling processes through 
randomly placed bridges is similar to the theory of metals 
with spin-flipping impurities [ p^ . In our case the role of 
the spin z component is played by the layer index. The 
main difference is that the different "spin" components 
are subjected to different (and uncorrelated) impurity 
potentials. We introduce the matrix notations for the 
diffusons D^^{Q,uj), where the indices k, 1, 111,71 —1,2 
indicate the layer index for each of the four incoming 
and outgoing legs. Such a matrix diffuson is a solution 
of the following Bette-Salpeter equation 
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where ({Q,uj) = + t^^ ~ {^Q^ ~ i^)^ ^'^id the 
combined tunneling and impurity scattering amplitude 
is given by 

W^n = -SklSmnSkm + — {l-5kl){l~ S,nn) ■ (9) 
T Tl2 

Since the elastic disorder in the two layers is uncorrelated 
one may look for a solution of Eq. in the form 

D'^niQ,^) = ASklS„,n5km + B{1 - &)(1 - (5„„0 , (10) 
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where U^^ is the matrix of bare interactions |6|-p0| and 
nfcn(Q, ^) is the polarization matrix which in the present 
case takes the form 



nfc„(Q,w) = 1/ [4„ + iwD^,'J(Q,w)] 



(14) 

Note, that by virtue of Eqs. ( p^ ) and (|Tl|) Ukn satisfy the 
particle conservation law: J2n ^kn{Q = 0,uj) — 0. 

Now, we are in the position to calculate Fxx{(^) for 
various values of uj. The second order contribution to 
the exchange transconductance (Fig. ^ comes from the 
expression (cf. Eq. ^)) Ui2D\\D\lDf^ = Ui2A^B. Ex- 
panding now up to the lowest order in hjTyi, one ob- 
tains Eq. (^. The fourth order originates from terms 
like UY\D\{D^^Df^ = UuAB'^, which should be taken 
in the lowest non- vanishing order in h/Ti2, and from the 
term Ui2DIIdIIdI'^ where Uu = -C/iini2C/22- As a 



result one obtains Fr 
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(w) with 
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The two terms on the r.h.s. of Eq. (|l|) are represented 
by the two classes of diagrams exemplified in Fig. ||. Sub- 
stitution of Eq. (15) in Eq. (^ and integration over fre- 
quency and slow momenta yields Eq. (3b). It is easy 
to see that this contribution is the dominant one at 

T < {Kd)h/Ti2. 

Finally in the limit of < ?i/ri2 one may neglect the 
second term in the expression for A and B [Eq. (^l])] and 
then Dll = DM = D\l = Df, = (1/2) {DQ^ ~ ilo)-\ 
This yields for uj < 1/ti2, Fxxi'-^) — FlfJ{Lo) with 



^^Q2^n+ 2^7,2 + ^22 



(16) 



then in the limit T12 3> r one finds 
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The function F^xi'^) which enters in expression (^) for 
the transconductance may be written as (cf. the dia- 
grams in Fig. [^) 



Fxxiu^) ^^Y^DQlY, UknDlt [D[^ + D[l] D, 

Q kiilrn 
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Here Ukn{Q,^) is the dynamically screened interaction 
between layers k and n. In the RPA approximation it is 
given by a solution of 



Eqs. (|l^ ) and (|lj) in the present case lead to U11 + U12 = 
U12 + U22 = iy^^{DQ^ - iuj)/DQ^. As a result of 
slow momenta summation one obtains for lo < h/Ti2: 



Fx'x^ (uj) = ie^ / {A-Kuu}). After substitution of Fx)c' {^) into 



Eq. (|6|) and integrating over the frequency one obtains a 
logarithm [Eq. (|3^)], similar to the logarithm in Ref. [p^ . 
For finite size systems the slow momentum integration 
should be replaced by a discrete summation. This re- 
sults in placing the low-frequency cutoff, lo ~ fiD/L'^, 
in Eqs. (|^), ( p^ ) and (16). As a consequence the tem- 
perature dependence of Eqs. (^-(|3^ is flattened out 
at r < hD/Lp'. In order to find the exact exchange 
contribution to transconductance in the transitions be- 
tween the various regimes (including the temperature in- 
dependent parts) one may directly substitute Eq. ( p^ ) 
into Eq. (0) and evaluate the resulting double integral. 
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Finally we briefly discuss the ballistic limit, Ti/t < T. 
For such a temperature range one should consider the di- 
agrams depicted in Fig. |l| without the diffusons. In this 



case the function F^x — Fxx 
4t2 



(bal) 



F^x\u^) 



Ui2{Q) 



is given by 
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where X = it^ + l/T), and C/i2(Q) is the screened 

interlayer interaction in the clean limit |]6|-|l0|. The slow 
momenta summation and frequency integration accord- 
ing to Eq. lead to a negative and temperature inde- 
pendent transconductance, which is given by: 



2 

One should not overestimate the significance of this neg- 
ative constant exchange contribution because of the ex- 
istence of two additional "parasitic" effects. Both con- 
tribute positive temperature independent constants to 
the measured transconductance. 

The first effect is related to a possibility of momentum 
conserving tunneling. In this case the transconductance 
is nonzero (and positive) at zero temperature even with- 
out any e-e interactions. It originates from the delocal- 
ization of the wavefunctions between the two layers, and 
may be described in the framework of the single-particle 
picture ]l7| , p^ . This type of (tunnehng) transconduc- 
tance may be differentiated from other mechanisms by 
a strong sensitivity to an in-plane magnetic field or a 
gate voltage. A Fermi-surface mismatch introduced by 
these factors leads to a rapid suppression of the transcon- 
ductance [TtIJI^ ]. The second effect originates from the 
classical (Hartree) interactions between electrons. In the 
presence of many bridges it may be visualized by a net- 
work of classical resistors. A simple calculation based 
on Kirchoff's laws leads to a positive transconductance. 
For L > \JDt\2 the network model predicts that ap- 
proximately half of the current applied to the first layer 
eventually leaks out to the second one. In this case the 
transresistance practically reduces to the resistance of a 
single layer of a doubled width [thus it is not surprising 
that we have recovered the Altshuler-Aronov interaction 
correction at low temperatures and at a large system size, 
Eq. (^)]. However, we stress that both mechanisms men- 
tioned above lead to a temperature independent effect, 
which is easily distinguishable from the strong temper- 
ature dependent exchange contribution to the transcon- 
ductance. 

To conclude we have examined the influence of point- 
like tunneling between the two layers, in combination 



1 VpT^ 



with e-e interactions, on the drag effect and the transcon- 
ductance. We found a distinguishable exchange contri- 
bution to the transconductance, which is negative, has 
a strong temperature dependence and is non-vanishing 
at zero temperature. The magnitude of the effect under 
optimal conditions may reach a fraction of /2tt% (for 
the transconductance). 
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